The synchronizability of highly clustered scale-free networks 
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In this letter, we consider the effect of clustering coefficient on the synchronizability of coupled 
oscillators located on scale-free networks. The analytic result for the value of clustering coefficient 
aiming at a highly clustered scale-free network model, the Holme-Kim is obtained, and the rela- 
tionship between network synchronizability and clustering coefficient is reported. The simulation 
results strongly suggest that the more clustered the network, the poorer its synchronizability. 

PACS numbers: 89.75,-k, 05.45.Xt 



Many social, biological, and communication systems 
can be properly described as complex networks with 
nodes representing individuals and edges mimicking the 
interactions among them 0,00. Examples are numer- 
ous: these include the Internet, the World Wide Web, so- 
cial networks, metabolic networks, food webs, and many 
others |j, 0, la, Q ■ Recent empirical studies indicate that 
the networks in various fields have some common char- 
acteristics, the most important of which are called small- 
world effect "s^ and scale- free property Networks of 
small- world effect have small average distance as random 
networks and large clustering coefhcient as regular ones. 
And the scale-free property means the degree distribu- 
tion of networks obeys the power-law form. 

One of the ultimate goals of researches on complex 
networks is to understand how the structure of complex 
networks affects the dynamical process taking place on 
them, such as traffic flow [Tolllll[T^IT3| . epidemic spread 
[iJilHillSi cascading behavior [13,13 13: and so on. In 
this letter, we concentrate on the synchronization, which 
is observed in a variety of natural, social, physical and 
biological systems [23, |2l|, 123 • The large networks of cou- 
pled dynamical systems that exhibit synchronized state 
are subjects of great interest. Previous studies mainly 
focus on the Wastts-Strogatz j^] networks and Barabasi- 
Albert networks, which have demonstrated that scale- 
free and small-world networks are much easier to syn- 
chronize than regular lattices ,23, 24,_25,,^]. Since many 
real-life networks are scale- free small-world networks, and 
there are already some models that can simultaneously 
reproduce the small-world and scale-free characteristics 
[27L 125 I29I l30l |. to investigate the synchronizabihty of 
the scale-free small-world networks if of great interest 
and importance. 

Since the scale-free networks are always of very small 
average distances [sj , the scale-fee small- world networks 
can also be referred as highly clustered scale-free net- 
works. One of the earliest highly clustered scale-free 
models is the Holme-Kim (HK) model [23|, which has 
successfully reproduced the indirectly acquainting mech- 



anism in real networks thus is closer to reality than BA 
model. In this letter, we will firstly give the analyti- 
cal result about the clustering coefficient of HK model, 
which is helpful for understanding the underlying evolu- 
tion mechanism of HK model. And then, we will investi- 
gate the relationship between network synchronizability 
and clustering coefficient based on HK model. 

As a remark, previous studies mainly concentrate 
on how the average distance and heterogeneity of de- 
gree/betweenness distribution affect the network syn- 
chronizability 2A^ ji2, M^, J^], while there are 
few systemic works about the effect of clustering coeffi- 
cient. Although there are a number of highly scale-free 
models, the HK model is a typical one who has tunable 
clustering coefficient thus provides us a good researching 
stage. This is the reason why we choose HK model as 
our theoretic template. 

The HK network is generated by the following pro- 
cesses: 

(1) In each step, m edges are added in the networks, 
and t is a discrete parameter which denotes the global 
time that the system totally goes. 

(2) First an edge is add with the probability H(fci), 



H(fc, 



(1) 



(3) Then, in the following to — I time steps, do a PA 
(preferential attachment) step with the probability p or 
a TF (triad formation) step with the probability 1 — p 
(see Ref. [27] for details). 

By using the rate-equation |37|. one can obtain the 
evolution of nodes' degree as follows: 



dh{t) 
dt 



= {l-(l- 



{l-(l- 



_fc^(^^_2)_ 



ki{t — l) -^mTF 1 ^2) 
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where Q,i denotes the set of neighbors of node i, ki{t) 
is the degree of node i at time step t. In the above 
formula, to^f is the number of edges that is connected 
following the rule of triad formation in each step while 



mpA denotes the number of those connected following 
the rule of preferential attachment. Denote ki ki{N), 
where is the network size, the two terms in the right 
side of the above formula are the degree increment rate 
of node i in PA and TF step, respectively. By using 
the initial condition ki(ti) — m, and the expressions of 
mpA = {m — 1)(1 —p) + l and uitf = pirn — 1), one can 
obtain the solution as follows: 



h{t) ^ ki{t;ti,m), (3) 

from which and by using of the continuum theory |36j| . 
one can get 
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(1) (2) (3) 



FIG. 1: The three parts of this figure represent three mech- 
anisms as: (1) node i is connected by a PA step while the 
neighbors of it is also connected by a TF step; (2) one of 
node i's neighbors is connected by a PA step while the node 
i is connected by a TF step; (3) node i as well as one of its 
neighbors is connected by a PA step respectively. 



P(fc,p) 



2m2 



+ A 



k^N' 



(4) 



where A is a quadric polynomial of p. Clearly, p{k) = k^^ 
in the limit case N oo. 

The clustering coefhcient of the whole network is the 
average of Ci over all nodes i, where Ci is the ratio between 
the number of edges among node I's neighbors which is 
denoted by rii and the total possible number. Then, 



2n,; 



ki{ki 1) 



(5) 



Using the rate-equation approach (37|] , the detailed ex- 
pression of c(fc), which denotes the average clustering co- 
efficient over all the fc-degree nodes, should be as follows: 
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Here we assume that toq (the number of initial nodes) is 
equal to m (the edges added each time step). The three 
items in the right side in the above expression is got from 
three mechanism shown in Fig. 1. 

Thus the clustering coefficient C can be solved as a 
function of the free parameter p. 



z{p) = 



c{k)P{k)dk. 



(7) 



In the above formula, kmax — > 2m\/N and kmin = 
|38|. For p G [0,1], c{p) can be simplified in linear ap- 
proximation. 



c{p) = B{Tn, N) + c{m, N)p. 



(8) 



FIG. 2: Clustering coefficient c(p) vs the parameter p in HK 
model. The square-dot curve represents the simulation result 
while the circle-dot curve denotes the analytic solution. The 
network size is A = 5000, and m = mo = 10. 



The extensive simulation results with different p and c for 
networks of different sizes strongly support the analytic 
results, especially in the larger-size networks, as shown 
in Fig. 2. 

Other simulations have been done about the variation 
of average path length I and standard deviation of degree 
distribution a. The results show that both the average 
path length I and standard deviation of degree distribu- 
tion a behave slightly variation with vary p. The intu- 
itionistic explanation is quite easily understood that a is 
directly related with degree distribution, that is to say, 
it is determined totally by the degree distribution, and 
hardly will the degree distribution vary when is N large 
enough. So we have the reason to neglect the effect of 
the two structural properties on synchronizability. 

Here, we concern the system of linear coupled limit- 
cycle oscillators on HK networks. Describing the state of 
the zth oscillator by Xi, the equations of motion governing 
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the dynamics of the N coupled oscillators are: 



N 



3 = i 



(9) 



where Xi = F{xi) characterizes the dynamics of individ- 
ual oscillators, G{xj) is the output function, K denotes 
the coupling strength and the N x N coupling matrix M 
is 



—ki, for i — j 
Mu^{ 1, forjsAi 
0, otherwise. 



(10) 



In the above expression, denotes the neighbors of node 
i. Because of the negative semidefinition and the zero 
sum of each raw of the matrix, all its eigenvalues are 
nonpositive real values and the biggest eigenvalue Aq is 
always zero. Thus the eigenvalues can be ranked as Aq > 
Ai > ... > \ni and Aq = Ai = if and only if the 
network is disconnected. 

In our coupled dynamic network, all the oscillators are 
identical and the same output function is used, the cou- 
pling fashion ensures the synchronization manifold is an 
invariant manifold and the nodes can be well approx- 
imated near the synchronous state by a linear opera- 
tor. Under these conditions, the eigenratio R = ^ can 
be used to measure the network synchronizability: the 
smaller it is, the stronger the synchronizability pSl Is^ . 

H IM IM EE El El El El El El 

We take only synchronizability R{p) and clustering co- 
efficient cij)) into consideration. Having simulated Rij)) 
for different configurations versus p (see Figure 3), we 
know that R{p) is positively correlated with p, and so is 
c(p) although it is not completely linear with p. Besides, 
we can easily find that when p is not very large, the curve 
is approximately linear (see Figure 4). 

Furthermore, we report the relationship between syn- 
chronizability R{p) and clustering coefficient c{p) as 
shown in Fig. 5. From the above figure, one important 
fact the curve reveals is exponential growing tendency. 
The function of fitted curve can be set as follows: 



i?(c) = A{m, N) + B{m, N)e~" 



(11) 



The terms A{m, N), B{m, N), and T{m,N) can be ob- 
tained by simulation. Here the simulation averaged over 
50 different realizations is to measure the effect of ran- 
dom fluctuation of degree distribution, since the degree 
distribution can both affect I and a. In our simulation, 
due to the average effect, the two parameters, I and a 
do not vary significantly so that we can only focus on 
clustering coefficient exclusively which causes the change 
of synchronizability. 

Compared to the previous works, the advantages of 
HK model is that the clustering coefficient can be tuned 




FIG. 3: Synchronizability R measured by -j^ vs the param- 
eter p in HK model. The simulation has been done under the 
condition that the size of the networks is ti = 1800. Other 
parameters are m = mo = 10. 
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FIG. 4: Clustering coefficient c{p) vs parameter p under the 
condition of A'^ = 1800, and m = mo = 10. The linear fitted 
line is for the first seven points. 



while the other structural properties are almost kept 
fixed. Therefore, combining the behaviors of c(p) ver- 
sus p and R{p) versus p, we obtain the relationship be- 
tween synchronizability and clustering coefficient. Fig. 
4 demonstrates that the larger the clustering coefficient, 
the poorer the synchronizability. Due to that the syn- 
chronizability R is determined by the ratio of maximal 
and minimal eigenvalues of coupling matrix which is ex- 
clusively related to the network topology, moreover c is 
the only varied topological property, thus clustering coef- 
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FIG. 5: Synchronizability R measured by ^ as a function 
of clustering coefficient c in the case of A'' = 1800, and m — 
mo = 10. The fitted curve follows the tendency of exponential 
growth. 



ficient plays a crucial role in synchronizability. The size 
effect has not been discussed so far in this letter. As 
the size N increases the average distance will get larger, 
as a result, the synchronization will become harder, but 
there will presumably not occur any qualitative changes, 
especially on the negative correlation between cluster- 
ing and synchronizability. Preliminary simulation results 
strengthen this conjecture. 

To ascertain the effect of each structural characteris- 
tics on synchronizability is a meaningful work because if 
future study can ascertain the relationships between each 
typical structural characteristics (Z, a etc.) and synchro- 
nizability exclusively, we can finally get the expression of 
synchronizability as function of those properties. Thus 
whether a network can achieve synchronization, and how 
is synchronizability of a specific structure can be easily 
predicted only by the topological characteristics of it. 
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